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ABSTRACT 


The  aerodynamic  loading  for  deformed  wings  with  elevons  in  both 
subsonic  and  supersonic  flow  is  considered.  The  solution  procedure  falls 
into  the  potential  flow  category  with  appropriate  restrictions.  For 
subsonic  flow,  a  lifting  surface  Kernel  function  formulation  is  used  in 
which  the  local  pressure  loading  for  both  wing  and  elevon  is  determined 
simultaneously  in  a  semi-closed  summation  manner.  Cases  under  study 
included  gaps  between  wing  and  elevon  in  addition  to  arbitrary  wing- 
elevon  deformations.  Results  for  all  cases  compared  very  well  with 
experimental  data.  Experimental  data  taken  in  a  low  speed  wind  tunnel  is 
also  presented  for  a  cropped  delta  wing  and  rectangular  elevon  in  which 
the  wing-elevon  gap  was  the  primary  test  variable.  For  supersonic  flow, 
3-D  supersonic  theory  forms  the  basis  for  the  solution  procedure. 
Deformations  are  accounted  for  with  the  use  of  doublet  paneling  added  to 
the  basic  3-D  solution.  Results  agree  very  well  with  existing  experi¬ 
mental  data.  The  gapped  elevon  and  thick  wing  trailing  edge  problem  is 
also  addressed  with  satisfactory  results. 
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I .  INTRODUCTION 


Subsonic  Flow 

It  has  been  observed  that  failures  associated  with  the  control 
surfaces  have  occurred  on  certain  missile  configurations.  It  is  believed 
that  these  failures  have  occurred  because  of  adverse  loading  on  the 
control  surface  due  to  large  deformations  in  the  wing  surfaces.  A  pre¬ 
liminary  study’'"  was  completed  which  attempted  to  identify  the  aerodynam¬ 
ic  changes  which  occurred  because  of  these  deformations.  The  present 
study  addresses  the  problem  in  more  detail  including  various  wing-elevon- 
gap  planform  configurations,  subsonic  and  supersonic  Mach  numbers,  moderate 
wing  and  elevon  deformations  of  an  arbitrary  nature,  and  blunt  trailing 
edges  for  the  wing. 

In  recent  years,  major  efforts  have  been  made  to  model  multi-element 
airfoil  configurations  (wing-elevon  combinations)  and  to  predict  resulting 
loads  and  moments.  Most  of  the  research  efforts  concerning  these  multi¬ 
element  lifting  surfaces  have  been  two-dimensional  analyses  as  typified 

2  3 

by  the  early  work  01  Glauert  ’  in  1924  and  1927  and  more  recently  by  the 
work  of  Halsey.  ’  However,  the  case  of  the  wing,  air  gap,  control 

6 

surface  combination  has  generally  been  neglected  as  pointed  out  by  Ashley 

who  also  expressed  the  need  for  such  an  analysis. 

Extension  of  analyses  to  include  three-dimensional  effects  or  finite 

wings  has  generally  employed  vortex  lattice  formulations  or  some  form  of 

7  8 

constant  pressure  paneling  such  as  that  of  Woodward  and  Lan.  .  The  Lan 
approach  has  received  considerable  attention  of  late  since  the  leading 


edge  suction  terms  and  Kutta  condition  at  the  trailing  edge  are  properly 

9 

accounted  for  in  the  analysis.  DeJarnette  extended  the  "strip"  approach 
of  Lan  to  provide  for  a  "continuous"  loading  in  the  spanwise  direction. 
Aside  from  the  Lan-De Jarnette  formulation,  general  vortex  lattice 
approaches  suffer  from  traditional  shortcomings  such  as  large  computer 
storage  and  long  run  times  even  on  today's  high-speed  computers.  Formu¬ 
lation  of  the  problem  using  the  Kernel  function  approach  in  its  usual  form 
as  typified  by  Cunningham^’ ^  has  been  generally  neglected  but  offers  the 

potential  for  increased  accuracy,  lower  computer  run  times  and  less 

8  9 

storage  than  does  the  vortex  lattice  approach.  The  Lan  -DeJarnette 
analysis  requires  comparable  run  times  and  storage  and  offers  a  viable 
alternative  to  the  solution  of  the  problem  at  least  for  single  element 
airfoils.  Extension  to  multi-element  airfoils,  as  far  as  is  known,  has 
not  been  completed. 

For  the  most  part,  previous  theoretical  analyses  have  been  restricted 
to  planar  lifting  surfaces  with  no  camber  or  twist.  Wing-elevon  con¬ 
figurations  general.lv  have  been  sealed-gap  cases.  However,  White  and 
12  13 

Landahl  ’  have  developed  a  procedure  requiring  the  method  of  matched 
asymptotic  expansions  which  is  used  to  determine  the  load  distribution 
when  a  gap  exists  between  the  wing  and  elevon  in  two-dimensional  flow. 

Both  the  vortex  lattice  and  Kernel  function  approaches  have  enjoyed 
certain  successes  but  certainly  the  vortex  lattice  has  been  more  exten¬ 
sively  used.  There  are  certain  problems  which,  though  amenable  to  vortex 
lattice  solutions,  are  more  easily  handled  by  the  present  lifting  surface 
techniques;  namely,  that  of  wing  irregularities  such  as  twist,  camber, 
and  arbitrary  deformations.  Certainly  camber  and  twist  have  been 


I 


addressed  In  previous  efforts,  but  analyses  of  moderate  "arbitrary" 
deformations  of  a  lifting  surface  have  not  been  found  in  the  literature. 


The  present  work  is  concerned  with  the  wing-elevon  problem  for 

configurations  with  moderate  wing  deformations  such  as  camber,  twist, 

and  general  deformations  due  to  high  wing  loading  or  distortion  of  the 

surfaces  due  to  aerodynamic  heating.  Lifting  surface  theory  is  employed 

14 

with  the  solution  following  the  procedure  established  by  Purvis  and 
Burkhalter,  et  al,^  The  loading  function  over  the  multiple  lifting 
surface  is  defined,  and  the  Kernel  function  is  then  integrated  over  the 
surfaces.  For  this  work,  the  gap  distances  are  considered  to  be  small 
enough  that  vortex  rollup  is  assumed  negligible. 

Comparison  is  made  with  experimental  data  on  several  general  con¬ 
figurations.  For  the  gap  case,  data  was  obtained  from  low  speed  wind 
tunnel  tests  with  gap  distance  as  a  primary  test  variable.  Data  for 
several  elevon  deflection  angles  and  hinge  line  locations  are  shown  for 
each  configuration  at  numerous  angles  of  attack.  Experimental  data  is 
lacking  for  wing  load  deformations  at  subsonic  speeds  but  comparisons 
are  made  with  wings  which  have  ridj.^te  camber  an!  twist. 

Supersonic  Flow 

For  supersonic  flow  there  are  several  approaches  which  address  the 
problem  of  predicting  loads  cn  a  planar  wing.  These  methods  range  from 
two-dimensional  shock  expansion  theory  to  three-dimensional  wing  theory 
to  vortex  lattice  theory.  Some  approaches  are  useful  for  variably  swept 
wings  while  others  are  only  applicable  to  wings  with  straignt  leading 
edges.  The  basic  elements  of  the  major  supersonic  flow  theories  will  be 
summarized  in  the  following  paragraphs. 


Shock-expansion  theory  is  two-dimens ionally  restricted  to  attached 
shocks  but  it  is  general  and  can  be  used  for  thick  and  thin  2-D  bodies 
alike.  However,  because  the  results  are  difficult  to  express  in  concise 
analytical  form,  shock-expansion  theory  is  mainly  used  for  obtaining 
numerical  solutions  for  diamond  shaped  airfoils.  If,  however,  one 
assumes  a  thin  wing  in  steady  flow  at  a  small  angle  of  attack;  i.e., 
small  perturbations;  then  the  approximate  relations  for  weak  shocks  and 
expansions  may  be  used.  This  thin  airfoil  theory  results  in  a  simple 
analytical  expression  for  calculating  the  pressure  coefficient. 

Although  these  two-dimensional  theories  are  valuable,  they  are  not 
sufficient  to  correctly  model  three-dimensional  flow.  Two  of  the  first 
three-dimensional  theories  were  point  source  theory  and  line  source 
theory.  These  theories  were  adapted  to  thin  supersonic  airfoils  with 
straight  leading  edges  and  various  types  of  symmetrical  crossections . 

In  each  theory  the  velocity  potential  for  the  wing  is  obtained  by 
analytically  integrating  over  the  forecone  from  the  point  in  question. 
Either  subsonic  or  supersonic  leading  edges  may  be  handled  in  this  manner 
as  well  as  delta  wings  with  supersonic  leading  and  trailing  edges.  These 
theories,  however,  do  not  allow  for  the  interaction  between  the  upper 
and  lower  surfaces  of  the  wing  which  would  result  from  subsonic  leading 
edges  as  well  as  from  subsonic  wing  tips. 

Thin  wings  with  straight  subsonic  leading  edges  may  be  handled  by 
conical  flow  theory.  As  the  flow  properties  are  constant  along  rays 
emanating  from  the  vertex  of  the  wing  the  flow  is  two-dimensional.  The 
Tschaplygin  transformation  is  used  to  transform  the  small  perturbation 
equation  in  polar  coordinates  into  the  two-dimensional  Laplace  equation. 


Complex  variable  techniques  developed  for  incompressible  flow  are  then 
used  to  obtain  an  analytical  solution  for  the  pressure  coefficients. 

This  technique  has  also  been  extended  to  wings  with  straignt  supersonic 
leading  edges  and  also  to  wings  with  straight  subsonic  and  supersonic 
leading  edges  combined. 

All  of  these  techniques  are  basically  inadequate  in  predicting 

loads  for  a  deformed  wing;  a  deformed  wing  in  this  case  meaning  moderate 

arbitrary  deformations  which  include  camber  and  twist.  Most  are  amenable 

to  general  planform  shapes  but  are  usually  restricted  to  straight  leading 

and  trailing  edges  and  do  not  include  fuselage  effects.  A  more  recent 
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attempt  to  include  curved  leading  edges  was  completed  by  Carlson  in 
which  a  vortex  lattice  scheme  was  used  to  model  the  wing  planform.  Al¬ 
though  resulting  computations  were  highly  oscillatory,  a  "smoothing" 
routine  introduced  by  Carlson  produced  acceptable  results  for  thin  planar 
wings.  Another  recent  formulation^  of  the  problem  utilizing  supersonic 
line  sources  has  also  produced  good  results  on  curved  planform  shapes 
and  does  not  require  any  post  data  manipulation. 

The  solution  technique  used  in  the  present  research  is  fundamentally 
based  on  three-dimensional  supersonic  flow  theory  as  outlined  in  Refs.  18, 

15  and  19,  but  is  modified  to  account  for  wing  deformations.  Potential 
flow  is  still  assumed  so  that  solutions  to  the  potential  equations  may 
be  added.  The  deformed  wing  is  "overlaid"  with  a  doublet  paneling  sheet 
whose  strength  everywhere  is  determined  through  matrix  inversion  tech¬ 
niques  and  the  doublet  solution  is  then  added  to  the  basic  three-dimensional 
supersonic  solution.  This  is  accomplished  by  subdividing  the  wing  planform 
into  very  small  panels  over  which  the  pressure  is  assumed  constant. 


Integration  is  completed  over  each  subpanel  in  closed  form  and  results 
may  then  be  written  in  summation  notation.  The  vorticity  paneling 
produces  perturbation  velocities  which  account  for  deformations  in  the  wing 
surface . 


II.  SUBSONIC  FLOW-THEORY 

Fundamental  Equations 

The  appropriate  equations  and  solution  procedure  for  the  Kernel 
function  formulation  used  in  the  subsonic  analysis  are  outlined  in  this 
section.  The  equations  are  well  known  (Refs.  20  and  21)  and  the  solution 
procedure  is  developed  in  detail  in  Ref.  14. 

Compressibility  is  included  in  the  potential  equation  which  is 
written  as 


<Kx 
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The  downwash  at  an  arbitrary  point  (x,y)  in  the  z=0  plane  due  to  an 

infinitesimal  area  (dx  dy  )  of  a  lifting  surface  is 
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where  £  and  n  are  dimensionless  chordwise  and  spanwise  variables,  respec¬ 
tively.  These  variables  are  made  non-dimensional  by  the  definitions: 


e  = 


x-x  (y) 
LE  ' 


c(y) 


n  = 

b/2 


(2) 


The  functional  form  of  the  pressure  loading  coefficient  is  assumed  as 


AC  (£.0>  = 
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n=0 


c(n) 
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m=0 


B  sin(2m+l) 6 
nm 
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where  n  =  cos  6  and  N  and  M  represent  an  arbitrary  number  of  chordwise 
and  spanwise  control  points,  respectively.  It  is  readily  seen  that  this 
form  produces  the  required  leading  edge  square  root  singularity,  satisfies 
the  Kutta  condition  at  the  trailing  edge,  and  has  slender  wing  behavior 
at  the  wing  tips.  A  discussion  of  the  logic  behind  this  assumed  form  of 
the  pressure  loading  is  presented  in  Reference  14. 

Over  a  sufficiently  small  subpanel  of  the  wing  planform,  AC^  is 
essentially  constant  and  may  be  taken  outside  the  integral  of  eq.  (1). 

The  resulting  expression  may  be  evaluated  in  closed  form  to  give 

ac  (C,n) 

Aw(x,y)  =  — B -  [K(x  ,y  )  -K(x„,y  )  -K(x  ,y  )  +K(x  ,y  )]  (4) 

8tt  22  *  1  12  11 

where,  as  an  example,  K(x^,y^)  is  evaluated  in  the  form 

(x-x  )  +  /(X-X  )2  +  &2(y-y  )2 

K(  y )  = - i - - - — 

1  (y-yx) 

+  3  loge[3(y-y1>  +  /(x-x^) 2  +  62(y-y1>a  ]  (5) 

The  other  K's  of  eq.  (4)  are  evaluated  by  making  the  appropriate  sub¬ 
stitution  of  the  subscripted  variables  in  eq.  (5).  Also  in  eq.  (4), 

ACp  (C,ri)  is  evaluated  at  the  centroid  of  the  subpanel.  The  wing  and 
subpanel  are  illustrated  in  Fig.  1. 

For  small  angles-of-attack,  the  boundary  condition  imposed  is  the 
requirement  of  no  flow  through  the  wing;  i.e., 

w(x,y)  +  sin  a(x,y)  =  0  (6) 
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The  total  downwash  at  any  control  point  due  to  the  entire  wing  is  found 

/ 

by  summing  the  effect  of  each  subpanel  of  the  surface  and  Eq.  16)  becomes 

V)  Aw(X',y  )  +  sin  cc(x.,y.)  =  0  (7) 

V  1  .1  1  J 


where  (x^,y.)  are  the  control  points,  and  are  located  according  to  a 
cosine  distribution1^ 


c(y .) 


Xi  =  XLE(yi)  +  2^  (1  •  C°S  N+I^  ’  *  ~ 


yj  =  2  (1  “ cos  ’  j  = 


(8) 


The  set  of  simultaneous,  linear  equations  for  the  NxM  unknown  loading 

coefficients,  B  's,  defined  by  substituting  eq.  (3)  into  eq.  (4)  and 
nm 

the  result  into  eq.  (7),  may  now  be  solved,  and  the  loading  at  any  point 
on  the  wing  is  defined  by  Eq .  (3) 


Loading  Function  for  Wing-Gap-Elevon  C onf i g urat ion 

The  above  analysis  is  for  a  wing  configuration  only.  When  a  control 
surface  (elevon)  is  added  to  the  wing,  the  load  distribution  will  of 
course  change.  i'urther  changes  in  this  distribution  occur  when  a  gap  is 
introduced  between  the  wing  and  eleven.  In  the  present  work,  it  is  assumed 
that  the  elevon  span  is  equal  to  the  wing  span,  and  the  elevon  hinge 
line  is  perpendicular  to  tiie  wing  centerline.  With  these  stipulations,  the 
solution  procedure  for  the  basic  wing  is  paralleled  except  for  the  loading 
function . 

First,  for  the  sealed -gap  case  with  a  deflected  elevon,  a  two-function 
formulation  using  eq .  (3)  to  describe  the  pressure  distribution  is  used. 

The  nondimensional  chordwise  variable,  F,  of  eq.  (3)  is  redefined  and  is 


considered  separately  for  the  wing  and  eleven.  For  the  wing  and  elevon, 
respectively,  it  becomes 


X-XLF,(y) 

r  = _ h. _ 

w  cw(y)  +  c-E(y) 


JF. 


X'XLEE(y) 


(9) 


This  geometry  is  shown  in  Fig.  2.  With  these  definitions,  it  is  easily 
verified  that  the  leading  edge  singularity  is  satisfied  on  both  the  wing 
and  elevon,  a:  d  the  Kutta  condition  is  met  at  the  elevon  trailing  edge. 
These  singularities  agree  with  usual  assumed  loading  function  with  con¬ 
trol  surfaces  pointed  out  by  Landahl  in  Reference  IT.  At  the  wing 
trailing  edge,  the  Kutta  condition  is  not  satisfied,  nor  is  the  usual 
singularity  obtained,  but  rather  a  finite  value.  This  loading  is  used  to 
allow  a  closed  form  solution  to  the  sectional  lift  integral.  Error  due 
to  this  assumption  is  considered  to  be  on  the  order  of  error  due  to  the 
assumption  of  an  infinite  loading,  while  the  advantages  of  a  closed  form 
solution  are  considerable.  This  (chordwise)  distribution  is  shown  sche¬ 
matically  in  Fig.  3a. 

When  a  gap  exists  between  the  wing  and  elevon,  the  loading  function 
over  the  wing  and  elevon  is  defined  in  order  to  create  the  effect  of  two 
wings.  The  leading  edge  singularity  and  Kutta  condition  must  be  satis¬ 
fied  on  each  surface.  Equation  (3)  produces  this  distribution  by 
redefining  the  chord-wise  variables  illustrated  in  Fig.  2  as 


^W  = 


x-x  (y) 

TP  v J  1 

L  w 


cw(y) 


(10) 
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for  the  wing  and  elevon,  respectively.  This  loading  is  shown  in  Fig.  3b. 
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Figure  4.  Control  Point  Location  and  Corresponding 
Matrix  for  Wing-Gap-Elevon  Configuration. 
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Separating  eq.  (12)  into  the  contributions  for  the  wing  and  elevon,  and 
using  the  definition  for  the  pressure  coefficient,  the  above  equation 
becomes 


b/2  x  _  x 

TEu  TE, 

-<./[/  Acpud\'  *  f 


4C  dxjdy 


(13) 
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It  is  convenient  to  define  the  sectional  lift  coefficient  as 
TE 


cc, 
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(14) 
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Considering  the  sectional  lift  for  the  entire  configuration  to  be 


«5(y  )  =  cc  (y  )  +  cc  (y  ) 

1  o  £  W  °  yE  0 

the  lift  becomes 
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(15) 
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where 
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Combining  eqs.  (9)  and  (10),  the  dimensionless  chordwise  variable  can 
be  defined  by 
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where 
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With  these  definitions,  the  sectional  lift  becomes 
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Using  eq .  (3)  for  the  pressure  distribution,  the  integration  of 
eq.  (19)  may  be  performed  as  in  Reference  15  to  give 
NW  M 

cc„(y  )  =  Y  y  I  B  sin(2m+l)0 
Z  o  L  L  "  r'm 
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where  represents  the  pressure  loading  coefficients  obtained  from 
eq.  (11)  corresponding  to  the  wing  loading,  and  B  represents  the 
coefficients  of  the  elevon  loading.  Also,  in  eq.  (20)  for  the  wing 
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Note  that  for  the  gap  case,  £TE  =  1,  and  eqs.  (21)  and  (22)  reduce  to 


eqs.  (23)  and  (24)  with  n  instead  of  p. 
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Returning  to  eq.  (16)  and  using  the  dimensionless  spanwise  variable. 


H,  of  eq.  (2),  the  total  lift  is 

1 

L  =  q  cc^(n)  dn 


(25) 
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With 


the  definition  of  the  lift  coefficient,  and  for  a  symmetric  wing, 


eq.  (25)  becomes 
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cc£(n)  dn 


(26) 


Similar  manipulation  of  the  force  equations  leads  to  the  total  moment 

about  the  root  chord  leading  edge  (of  the  wing)  as 
1 


M 


-  / 

Sc  J 


2~  (n)  dn 


C  C 


(27) 


with  the  sectional  pitching  moment  coefficient  given  by 
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where  the  first  two  terms  are  seen  to  be  the  moment  due  to  the  wing,  and 
the  last  two  are  the  moment  due  to  the  elevon. 
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For  the  wing  alone,  the  integration  of  the  lift  and  moment  equations 
is  performed  analytically  in  closed  form.  However,  because  of  the 
discontinuity  resulting  from  the  addition  of  an  elevon,  eqs.  (26)  and 
(27)  must  be  integrated  numerically  when  a  gap  exists  or  the  elevon  is 
deflected . 

In  defining  the  lift  and  moment  coefficients  for  the  entire  configur¬ 
ation,  the  reference  area,  S,  is  the  total  surface  area  of  the  wing  and 
elevon  combined,  excluding  the  gap  area.  Similarly,  the  reference  length 
in  the  moment  coefficient  equation  is  the  mean  aerodynamic  chord,  c,  of 
the  wing  and  elevon,  without  considering  gap  distance. 

Elevon  Hinge  Moments 

When  an  elevon  is  included  in  the  configuration,  elevon  hinge  moments 

are  important  in  determining  control  forces  for  the  vehicle.  For  a 

straight  elevon  leading  edge  with  the  hinge  line  located  at  the  leading 

edge,  the  total  elevon  hinge  moment  is 
1 
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where  the  sectional  hinge  moment  coefficient  is  defined  as 
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Writing  eq.  (29)  in  terms  of  the  moment  coefficient  and  for  a  symmetric 
wing,  it  becomes 
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Considering  both  the  wing  and  elevon  contributions  of  eq.  (27)  separately, 
it  is  noted  that  eq.  (30)  is  the  elevon  moment  about  the  elevon  leading 
edge,  the  solution  of  which  is  given  by  the  last  term  of  eq.  (28),  or 
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sin(2q+l) 6 
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Substituting  eq.  (32)  into  eq.  (31),  the  integration  may  be  evaluated  in 
closed  form  (Reference  14)  to  obtain 
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(34) 


For  the  elevon  hinge  moments,  the  reference  area,  S^,  is  the  total 

area  of  the  elevon.  Also,  the  reference  length,  c  ,  is  the  mean  aero- 

E 

dynamic  chord  of  the  elevon. 

Thus  far,  the  analysi  has  be  n  limited  to  a  hinge  line  at  the  elevon 
leading  edge.  If  the  elevon  lift  is  also  .mown,  the  lift  *  id  moment 
produce  a  couple  at  the  elevon  leading  edge.  The  moment  about  any  chord- 
wise  station  may  be  found  by  summing  moments  at  that  point.  This  gives 
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where  x  is  the  chordwise  location  of  the  hinge  line  measured  from  the 
HL 


elevon  leading  edge.  Dividing  by  q  Sc  gives  the  equation  in  coefficient 

E  E 


form 
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where  C„  is  given  by  eqs.  (33)  and  (34). 

n 


The  total  lift  for  the  elevon  is 
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(37) 
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with  the  sectional  lift  of  the  elevon  given  by 
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c\<n)  ’  CE 


(38) 


It  is  readily  seen  that  this  is  the  second  integral  in  eq.  (19)  and 


the  solution  is  given  by  the  last  term  of  eq .  (20);  thus 
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where  I  is  given  by  eqs.  (23)  and  (24). 
P 


For  the  elevon,  the  lift  may  be  obtained  in  closed  form.  Substi¬ 
tuting  eq.  (39)  into  eq.  (37),  and  integrating  gives 
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where,  in  the  loading  coefficients  B  ,  °nly  the  q=0  term  appears. 


Thus,  using  eqs.  (40)  and  (33)  in  eq .  (36),  the  elevon  hinge  moment 
about  the  hinge  line  is  obtained  in  closed  form.  The  only  requirement 


irfMj-W? 
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for  the  analysis  is  that  the  hinge  line  be  perpendicular  to  the  wing 
root  chord. 

Non-Planar  Wings 

Lifting  surfaces  which  do  not  lie  entirely  in  the  x-y  plane  present 
a  different  view  to  the  flow  field  than  the  planar  wing  previously  con¬ 
sidered.  To  analyze  non-planar  wings,  it  is  necessary  to  consider  the 
angle  that  each  control  point  and  each  grid  element  makes  with  the  x-y 
plane.  Non-planar  effects  of  primary  interest  include  wing  twist,  cam¬ 
bered  airfoils,  and  arbitrary  ehordwise  and  spanwise  deformations. 

For  the  case  of  wi  lg  twist,  at  any  given  spanwise  location  on  the 
wing,  the  chord  line  lies  in  a  plane  that  differs  from  the  initial  x-y 
plane  by  the  twist  angle,  a,^,  at  that  station.  This  angle  may  be  in¬ 
cluded  in  the  boundary  condition  such  that  eq.  (6)  becomes: 


w(x.,y.)  +  sin  [a(x  ,y.)  +  a  (x  ,y.)]  =  0  (41) 

i  J  i  J  TW  i  j 

The  boundary  condition  for  the  eleven,  eq.  (11),  may  be  rewritten  in  a 
similar  manner. 

For  cambered  or  arbitrarily  deformed  airfoils,  the  chord  line  may 
be  replaced  by  the  mean  cam!  ev  line.  The  downwash  produced  by  each  grid 
element  is  proportional  to  the  angle  which  that  element  makes  with  the 
x-axis ,  or 


Aw(x.,y.)  =  Aw(x . ,y . ) cosa . 
1  J  1  J  1 


(42) 


where  is  the  angle  relative  to  the  x-axis  and  is  given  by 


ol  =  tan-1 
i  dx 


lx. 

1 


(43) 
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This  is  Che  velocity  component  in  the  z-f‘irection;  however,  each  control 
point  also  has  a  z-coordinat e .  The  boundary  condition  is  for  no  flow¬ 
through  the  wing,  or  that  the  normal  component  of  velocity  to  the  wing 
be  zero.  The  angle  that  the  control  point  makes  with  the  x-axis  is 
given  by 

Ot  =  tan- 1  (4^)1  (44) 

c  n  :  |x 

c 

Thus  the  component  of  velocity  normal  to  the  wing  is 

Aw(x.,y.)  =  Aw(x. ,y . )cosu  (45) 

1  j  1  j  c 

It  should  be  noted  that  the  boundary  condition  of  the  deformed  wing  is 
satisfied  on  the  x-axis,  although  the  coordinate  of  the  mean  line  lies 
off  the  axis;  however,  this  is  equivalent  to  the  planar  wing  approxi¬ 
mation  essential  to  thin  airfoil  theory.  The  actual  normalwash  is  found 
by  combining  equations  (42)  and  (45) 


Aw(x.,y.)  =  Aw(x.  ,y .  )cos( '..-a  ) 
1  j  i  j  ic 


( 46 ) 


By  defining  a  function  z(x)  many  variations  and  amplitudes  of  camber 
and  wing  deformation  in  the  chordwise  direction  may  be  obtained.  Com¬ 
binations  of  twist  and  camber  or  deformation  may  also  be  analyzed.  How¬ 
ever,  the  condition  of  no  flow  separation  still  applies  at  all  points. 
Interference  Effects 

In  the  present  analysis,  interference  effects,  such  as  fuselage  per¬ 
turbations  and  leading  edge  suction,  are  also  included.  An  infinite 
line  doublet  is  used  to  simulate  fuselage  interference  effects  on  the 
induced  velocity  over  the  wing.  The  potential  equation  for  the  line 


! 

1 


_  4 


doublet  may  be  differentiated  to  determine  the  velocity  component  normal 
to  the  plane  of  the  wing.  To  calculate  the  lift  and  moment  of  a  fuselage 
in  the  presence  of  a  wing,  an  image  wing  inside  the  fuselage  is  emploved. 
The  sectional  characteristics  of  the  actual  wing  are  integrated  over  the 
image  wing  span  to  determine  the  lift  and  moment.  Both  of  these  principles 
are  well  known  and  the  theoretical  development  can  be  found  in  Reference 
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Leading  edge  suction,  which  makes  a  significant  contribution  to  the 
lift  of  highly  swept  wings  at  large  angles  of  attack,  has  also  been 
incorporated  into  the  present  theory.  Purvis  (Ref.  23)  has  shown  that, 
for  an  assumed  pressure  distribution  in  the  form  of  eq.  (3),  and  for  an 
elliptic  spanwise  lift  distribution,  a  general,  closed  form  solution  for 
the  leading  edge  suction  is 

C2 


VLE 


=  (CL  sina 


cos') 


cost 

cosA 


where  C  represents  the  potential  flow 
LP 

by  eq.  (26).  Also,  from  Reference  23, 


lift  coefficient 
the  moment  may  be 


(47) 

which  is  given 
written  as 
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[  n/l-r,‘  +  sin  '  ( n)  ] dy 
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Using  the  definition  of  q  in  eq .  (2)  gives 


(48) 
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(49) 


For  a  constant  leading  edge  sweep  angle,  this  integral  may  be  evaluated 
in  closed  form.  The  coordinate  of  the  leading  edge  is  given  by 


=  r,  cosA 


therefore,  eq  .  (49)  becomes 

1. 106036C, 

tanA 


?'VlE  c 

Integration  produces 
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(51) 


Since  the  current  gap  analysis  assumes  small  gaps  such  that  there  is 
no  vortex  rollup,  interference  effects  will  be  confined  to  the  fuselage 
and  leading  edge  suction.  These  correction  terms  are  necessary  to  allow 
a  comparison  between  theory  and  experiment  in  many  instances;  therefore, 
by  including  them,  a  more  versatile  analysis  is  obtained. 

Order  of  Polynomials 

Before  proceeding,  comments  on  the  effect  of  polynomial  order  on 


the  solution  is  in  order.  The  unknown  constants,  B  in  eq .  (3)  are 

nm 


really  coefficients  for  the  various  terms  in  the  polynomial  describing  the 
pressure  load  on  tb-»  wing.  It  has  been  suggested  that  if  th.  polynomial 
order  is  increased  to  a  large  number  of  terms  (i.e.,  >  50)  that  the  sec¬ 
tional  coefficients  such  as  cc  will  tend  to  "oscillate"  in  the  spanvise 
direction.  For  the  formulation  used  in  the  present  analysis,  no  oscil¬ 
lations  are  observable  at  least  up  to  50  terms  in  the  equation.  However, 
oscillations  do  indeed  occur  if  terms  on  the  order  of  80  or  more  are 
used  in  the  pressure  loading.  It  does  not  appear  that  this  is  due  to 
"numerical"  instability  as  would  be  the  case  if  numerical  integration 


were  used,  but  rather  due  to  overkill.  That  is,  the  pressure  loading 
function  is  rather  general  and  for  most  wings  only  a  few  terms  (•'  20)  are 
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required  for  a  solution  If  more  terms  are  used,  it  is  analogous  to 
fitting  a  50th  degree  polynomial  to  a  linear  curve;  oscillations  are  the 


inevitable  results. 


III.  EXPERIMENTAL  MODEL  AND  TESTS  -  SUBSONIC 


Model 

To  verify  the  gap  effects  and  subsequent  elevon  loading,  an  exper¬ 
imental  model  shown  in  Fig.  5  was  fabricated  and  tested  in  a  low  speed 
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wind  tunnel.  The  baseline  configuration  consisted  of  a  thin,  low- 
aspect  ratio,  cropped  delta  wing  with  rectangular  elevons.  These 
surfaces  were  symmetrical  about  both  the  chord  line  and  the  fuselage 
centerline.  The  fuselage  was  included  to  facilitate  mounting  of  the 
hinge  moment  balance  mechanism.  An  ogive  nosecone  and  tailcone  were 
also  included.  The  fuselage  and  wing  tips  were  designed  to  allow  for  a 
two-inch  translation  in  the  elevon  mounting  position,  which  could  be  used 
to  vary  the  gap  distance  and/or  the  elevon  hinge  line  location.  The 
model  was  floor-mounted  to  an  external  six-component  pyramidal  balance. 
Test  Procedures 

The  model  was  tested  at  a  dynamic  pressure  of  3.5  inches  of  water 
corresponding  to  a  speed  of  approximately  126  feet  per  second  and  a 
Reynolds  number  of  0.76  million  per  foot. 

Data  was  obtained  for  an  angle  of  attack  range  from  -19  to  +5 
degrees  in  2-desree  increments.  The  negative  range  was  used  to  minimize 
strut  interference  effects  on  the  elevons.  Several  elevon  deflection 
angles  (0°  to  20°),  gap  distances  (0  to  75%  Cg)  and  hinge  line  locations 
(0  to  50%  cj  were  tested. 

ti 
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Figure  5.  concluded 


IV.  RESULTS  -  SUBSONIC  FLOW 


Results  for  a  wing  configuration  are  compared  to  other  well  known 
theories  in  Figs.  6  and  7.  Agreement  is  good  for  all  cases  considered. 

In  obtaining  results  using  the  present  theory,  the  wing  subpanel  mesh  of 
constant  loading  was  45  by  45  (chovdwise  by  spanwise) ,  and  the  control 
point  grid  was  4  by  4.  The  total  CPU  time  required  was  45  seconds  on  an 
IBM  3031  compared  to  1  minute  for  the  method  of  Lan  on  the  Honeywell 

Q 

635.  Storage  requirements  for  the  present  method  are  negligible  on  a 
large  machine.  This  mesh  size  may  be  reduced  such  that  total  CPU  time 
is  reduced  to  20  seconds,  with  a  change  in  accuracy  of  approximately  5 
percent . 

A  comparison  of  theory  and  experiment  for  the  "base"  configuration  of 
the  wing-fuselage  model  is  presented  in  Fig.  8.  Both  linear  theory  and 
leading  edge  suction  are  included.  The  importance  of  the  nonlinearity 
introduced  by  leading  edge  suction  in  the  high  angle  of  attack  region  is 
easily  realized  rrom  this  figure. 

The  theoretical  pressure  distribution  for  a  deflected  elevon  is  com¬ 
pared  with  experimental  data  (from  Ref.  25)  in  Fig.  9.  From  this  figure, 
it  is  seen  that  the  assumed  functional  form  for  the  chordwise  pressme 
distribution  for  the  sealed-gap  case  represents  an  adequate  model  of  the 
actual  flow.  Total  lift  and  hinge  moments  are  presented  in  Fig.  10 
where  agreement  is  good,  particularly  in  the  low  angle  of  attack  regime. 

The  analysis  developed  for  the  general  case  of  deformed  wings, 
elevons,  gap,  and  fuselage  is  used  to  predict  the  lift,  pitching  moment. 


Lan's  theory  (Ref.  8) 


Present  theory 
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Comparison  of  Two  Theories  for 
Swept  and  Tapered  Wing  of  Aspec 
Ratio  2.828. 
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and  hinge  moment  slopes  for  the  model  tested  in  the  low  speed  tunnel  and 
is  shown  as  a  function  of  gap  distance  in  Fig.  11. 

Certainly  the  most  interesting  trend  with  respect  to  the  gap  is  that, 
for  all  cases,  hinge  moments  increase  slightly  as  gap  distance  increases. 

It  is  also  noted  from  Fig.  11  that  the  coefficient  slopes  vary  more 
noticeably  with  changes  in  elevon  deflection  while  the  changes  with  angle 
of  attack  are  essentially  negligible.  Figure  12  illustrates  the  theoret¬ 
ical  spanwise  pressure  distribution  over  the  wing  and  elevon  for  the 
present  cropped  delta  wing  configuration.  This  figure  confirms  the 
predicted  increase  in  the  loading  over  the  elevon,  and  shows  the  slight 

loading  change  with  respect  to  angle  of  attack.  This  loading  form  is  in 

12 

agreement  with  the  two-dimensional  results  obtained  by  White,  ett  al,  as 
well  as  the  present  experimental  data.  In  Fig.  13,  the  hinge  moments  are 
presented  for  several  hinge  line  locations  on  the  elevon.  It  is  seen 
here  that  there  is  a  chordwise  point,  in  this  case  between  15  and  20  per¬ 
cent  of  the  elevon  mean  aerodynamic  chord,  where  the  hinge  moments  are 
independent  of  gap  distance.  It  is  noted  that  this  point  does  not  coin¬ 
cide  with  C  =0,  which  gi  '^s  the  elevon  a<-r^dynam/ c  center.  This  aspect 
H<5 

could  be  significant  in  the  design  of  vari'ble  geometry  wings  where  gap 
distances  are  subject  to  change.  Also  the  minimum  value  of  the  hinge 
moment,  which  corresponds  to  minimum  control  forces,  occurs  for  the  sealed- 
gap  case. 

The  slopes  of  the  pitching  moments  curves  illustrated  in  Fig.  11  are 
seen  to  be  essentially  constant  in  nature.  Since  the  moment  is  referenced 
to  a  fixed  point  on  the  wing,  the  increase  in  the  moment  arm  with  increas¬ 
ing  gap  coupled  with  increased  elevon  loading  might  tend  to  suggest  an 
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Figure  11.  Aerodynamic  Characteristics  for  a 
Gapped  Configuration. 


increase  in  the  n.oment.  However,  from  Fig.  12  it  can  be  seen  that  the 
change  in  pressure  distribution  more  than  offsets  these  factors.  Yet, 
as  the  gap  becomes  large,  the  theoretical  pressure  distribution  must 
approach  that  of  two  separate  wings,  without  interference.  The  moment 
referred  to  a  fixed  point  on  the  wing  must  become  infinitely  large  as 
the  gap  distance  approaches  infinity. 

For  large  gaps,  it  was  found  experimentally  (Figs.  14  and  15)  that 
there  is  a  considerable  increase  in  the  total  lift  and  moment  for  an 
elevon  deflection  of  -20°.  The  decrease  in  the  forces  becomes  larger 
with  increasing  gap  distance.  However,  elevon  hinge  moments  did  not 
reflect  this  trend,  indicating  that  it  is  due  more  to  the  gap  rather  than 
to  flow  separation  over  the  elevon.  In  any  case,  it  is  apparently  a 
viscous  phenomena  and  cannot  be  modeled  properly  using  potential  flow 
techniques . 

The  experimental  data  for  different  hinge  lines  is  presented  in 
Figs.  16  and  17.  Theoretically,  there  are  no  changes  in  the  total  forces 
and  moments  with  a  change  in  the  hinge  line.  This  result  also  appeared 
in  the  experimental  data  despite  the  fact  that,  for  a  deflected  elevon,  the 
elevon  leading  edge  no  longer  lies  in  the  plane  of  the  wing.  From  the 
data  in  Fig.  16,  it  is  seen  that  the  hinge  moments  aL  x^  /c  =  .254,  or 
approximately  the  quarter  chord  of  the  elevon,  are  quite  small  and  show 
little  variation  with  angle-of-at tack  or  elevon  deflection.  This  in¬ 
dicates  that  this  point  is  near  the  aerodynamic  center  for  the  elevon. 
However,  the  slopes  change  with  a  change  in  gap  distance,  thus  the  aero¬ 
dynamic  center  must  also  be  gap  dependent,  as  previously  shown  on  Fig.  13. 
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Results  for  arbitrarily  deformed  airfoils  are  difficult  to  find  in 
the  literature;  therefore,  perhaps  the  best  basis  for  comparison  is  two- 
dimensional  thin  airfoil  theorv.  Table  1  illustrates  theoretical  results 
for  the  two  arbitrary  nean  lines  shown.  The  present  theory,  extended 
to  two-dimensional  results  and  thin  airfoil  theory  show  good  agreement. 

For  cambered  airfoils,  experimental  data  may  be  obtained  from  Ref.  26. 
Results  are  presented  for  the  NACA  747A315  laminar  flow  airfoil  in  Fig.  18. 
The  current  theory  predicts  experimental  data  very  well.  Again,  the  thin 
airfoil  and  present  theory  are  closely  matched.  This  result  should  be 
expected  since,  foe  two-dimensional  flow,  the  present  theory  would  become 
thin  airfoil  theory  as  the  limit  of  the  wing  subpanel  areas  with  a  constant 
ACp  approached  zero. 

For  a  finite  wing.  Fig.  19  compares  theoretical  predictions  with 
experimental  data  (Ref.  27)  at  several  Mach  numbers  for  a  wing  that  is  both 
cambered  (V  C  \  a  -1.0  mear  <  .mbei  !  !ne)  and  <-w’'sted.  It  should  be  noted 
that  the  e.tper-f  mcr  tal  model  was  tapered  and  t!e  maximum  camber  varied  with 
the  spanwise  oorciinate.  In  the  theoretical  model  the  wing  was  replaced 
by  the  mean  camber  line  at  each  20  percent  of  the  span  and  this  distribu¬ 
tion  remained  constant  for  that  section.  Also,  a  linear  spanwise 
distribution  of  twist  was  assumed,  while  there  was  a  slight  variation 


from  this  experimentally.  Nontheless,  agreement  of  theory  and  experiment 
is  very  good  over  the  entire  experimental  Mach  number  range.. 


Table  1.  Results  for  Arbitrary  Mean  Tines. 


V.  CONCLUSIONS  -  SUBSONIC  FLOW 

From  the  preceeding  results,  it  is  seen  that  the  present  theory 
produces  good  results  in  the  inviscid  flow  regime.  Configurations 
ranging  from  a  basic  wing  to  varying  combinations  of  wings,  elevons,  gap, 
and  fuselage  may  be  modeled  as  well  as  moderate  wing  and/or  elevon  defor¬ 
mations.  Since  compressibility  effects  are  included,  reasonable  results 
up  to  large  (subsonic)  Mach  numbers  are  generally  obtained. 

For  the  finite  gap  case,  the  theory  provides  an  excellent  means  for 
obtaining  early  design  characteristics  in  a  short  time  for  small  angles- 
of-attack  and  elevon  deflections.  However,  it  is  shown  that  for  large 
elevon  deflections  there  was  a  considerable  decrease  in  total  lift  and 
moment  characteristics  when  a  large  gap  was  present.  This  appears  to  be 
a  viscous  phenomena  and  therefore  cannot  be  modeled  with  potential  flow 
techniques.  Vortex  rollup,  which  must  be  included  for  very  large  gaps, 
has  also  been  neglected  in  the  present  theory.  Results  for  cambered, 
twisted,  and  arbitrarily  deformed  wings  may  also  be  obtained  easily  with 
the  present  theory.  Results  are  good  as  long  as  the  requirements  of 
potential  flow  are  met  on  the  surface. 

Several  advantages  of  the  current  method  occur  as  a  consequence  of 
the  assumption  of  a  constant  AC^  over  a  sufficiently  small  subpanel  of  the 
wing.  The  wing  can  be  divided  into  many  subsections,  yet  only  a  relative¬ 
ly  small  matrix  must  be  inverted  to  obtain  the  load  distribution.  This 
assumption  also  allows  AC^  to  be  taken  outside  the  downwash  integral  and 
yields  a  closed  form  solution  to  the  equation.  As  a  result,  this  method 
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VI.  SUPERSONIC  FLOW  -  THEORY 


Fundamental  Equations 

The  theoretical  approach  for  the  supersonic  Mach  number  regime  is  a 
potential  flow  formulation  in  which  solutions  are  superimposed.  The 
source  solution  introduced  by  Eward^S  forms  the  basis  of  3-D  supersonic 
theory  used  in  the  present  analysis  and  vorticity  paneling  is  added  to 
account  for  wing  deformations.  The  vorticity  paneling,  as  previously 
introduced,  is  actually  a  continuously  distributed  doublet  panel  which 
satisfies  the  supersonic  potential  equation 
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As  for  subsonic  flow,  the  vorticity  distribution  y(x  ,y  )  may  be  replaced 

o  o 

by  the  pressure  loading  coefficient  JjAC  (x  ,y  ).  The  downwash  equation 

poo 

is  obtained  by  differentiation  in  the  standard  manner  and  when  evaluated 
in  the  z=0  plane  (planar  flow)  becomes 
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where  the  area  Sq  is  that  contained  on  the  wing  surface  in  the  upstream 
running  Mach  cone  emanating  from  the  field  point  (x,y).  Although  the 
vorticity  paneling  is  formulated  for  planar  flow,  it  may  be  used  to 
account  for  non-planar  velocity  components  if  the  slope  of  the  nonplanar 
surface  is  properly  taken  into  account.  A  more  detailed  discussion  of 
how  this  is  done  is  explained  later. 
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The  assumed  pressure  loading  consists  of  the  basic  three-dimensional 
supersonic  flow  term  plus  the  perturbation  term  from  the  vorticity  panel¬ 
ing  and  is  then  written  as 


N  M 

AC  (x  ,y  )  =  P(x  ,y  )sin(al0  )  +  P(x  ,y  )  £  J  B. 
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The  first  expression  on  the  right-hand  side  is  the  3-D  loading  term 
modified  to  account  for  local  wing  deformations  and  is  valid  for  both 
subsonic  and  supersonic  leading  edges.  The  second  term  uses  the  3-D 
pressure  functions  PCx^.y^)  as  a  weighting  function  which,  when  multi¬ 
plied  by  the  unknown  coefficient  B  ,  and  the  series  p1”  cos  accounts 

Lt  Z. 

for  small  deformations  in  the  wing  surface  in  the  upstream  Mach  cone  from 
the  point  (x,y).  The  order  of  the  terms  in  the  second  expression  seldom 
exceeds  3  or  4  and  for  planar  wings,  the  entire  second  term  is  unnecessary. 
Increasing  the  polynomial  order  for  larger  wing  deformations  or  higher 
angles  of  attack  does  not  seem  to  provide  better  agreement  with  experi¬ 
mental  data. 

The  form  of  the  assumed  pressure  leading  term  to  account  for  per¬ 
turbation  velocities  was  chosen  as  a  polynomial  in  the  spanwise  direction 
and  a  trigometric  or  Fourier  series  type  function  in  the  chordwise 
direction.  No  strong  justification  for  these  assumptions  can  be  given 
except  that  in  the  spanwise  direction  the  3-D  loading  term  should  produce 
acceptable  results,  even  for  deformed  wings,  and  corrections,  which  the 
polynomial  provides,  should  be  small.  Subsequent  good  agreement  with 
experimental  data  seem  to  justify  this  assumption.  In  the  chordwise 
direction  where  more  severe  changes  occur  in  the  pressure  loading,  a 
trigometric  correction,  which  is  somewhat  more  general  than  the 


polynomial,  was  chosen  with  subsequent  good  agreement  with  experimental 
data. 


As  in  subsonic  flow,  the  ACp  term  is  removed  from  the  integral  in 
eq.  (53)  and  the  resultant  expression  is  integrated  in  closed  form  over 
a  small  subpanel.  Results  of  this  integration  over  various  shaped 
subelements  is  presented  in  Ref.  15. 

It  is  convenient  to  define  the  value  of  this  integral  as  BAK  so  that 
the  total  downwash  at  a  point  (x,y)  becomes 


w(x,y)  =  l  (AK).[AC  (xQ,yo) 
i=l 


I.  (f-) 

1  4tt 


(55) 


where  I  is  the  total  number  of  subpanels  in  the  region  (cone)  of  integra¬ 
tion.  The  boundary  condition  is 

w(x,y)  +  sin  (a+0^)  =  0  (56) 

where  0  is  the  local  deformation  angle  at  a  control  point.  Rearranging 
c 

and  combining  equations  (54) ,  (55)  and  (56)  yields 
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(57) 


where  L  =  n  .+  mN. 

Body  perturbation  effects  on  the  wing  loading  are  accounted  for  by 
assuming  that  the  body  is  cylindrically  shaped.  The  cylinder  is  gen¬ 
erated  with  an  infinite  line  doublet  whose  axis  is  in  the  z-direction 


perpendicular  to  the  wing  planform.  The  results  of  adding  the  body  to  the 


analysis  simpiy  adds  another  term  to  the  right-hand  side  of  eq. 
after  a  little  manipulation,  becomes 
4na2  sina 

6<a+n)2 


(57)  which 


(58) 


where  a  is  the  nondimensional  body  radius  defined  by  a  = 
Equations  (57)  and  (58)  are  combined  and  rearranged 


rfa/(b/2). 

as 
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(59) 


The  above  equation  represents  a  set  of  N  simultareous  linear  alge¬ 
braic  equations  and  is  conveniently  arranged  in  matrix  form  for  a 
Gaussian  reduction  solution  for  the  unknown  coefficients  B  .  Once  the 

L 

unknown  coefficients  are  determined,  the  localized  loading  is  obtained 
from  equation  (54) . 

Lift  and  Pitching  Moment  for  Wing  Elevon  Combinations 

In  the  preceding  section,  the  basic  equations  were  presented  for 
finding  the  resulting  pressure  distribution  on  a  wing  or  elevon  in  super¬ 
sonic  flow.  For  the  wing  alone,  these  equations  may  be  applied  directly 
since  there  is  no  upstream  influence  due  to  the  elevon.  After  the 
pressure  loading  coefficients  in  eq.  (54)  are  determined,  the  pressure 


coefficient  may  be  computed  at  any  point  on  the  wing  surface  and  appro¬ 
priate  subsequent  integration  produces  the  lift. 
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For  the  elevon,  a  similar  procedure  is  followed  with  the  exception 
that  the  elevon  is  flying  in  the  wake  of  the  wing  and  therefore  the  down- 
wash  produced  by  the  wing  and  the  elevon  must  be  considered.  The  mechanism 
used  to  compute  this  downwash  is  discussed  in  another  section. 

The  lift  for  the  wing-elevon  configuration  is  the  sum  of  that  pro¬ 
duced  by  the  wing  and  that  produced  by  the  elevon.  It  is  convenient  to 
define  nondimens ional  chordwise  and  spanwise  variables  for  the  wing  as 


= 


(x-x  )/C 
LE  w 


(60) 


and 

n  =  y/(b/2)  (61) 

W 

For  the  elevon,  using  the  same  coordinate  system,  these  definitions 
become 


CE  =  (x-CR-e)/CE  (62) 

and 

nE  =  y/(b/2)  (63) 


With  these  definitions,  integration  of  the  pressure  distribution 
over  the  wing  for  the  lift  becomes 


ct  =  ~  /  f  (Ac  c  dc  )dn 

L  2S  /  /  Pwww 
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and  for  the  elevon 
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It  is  convenient  to  define  sectional  properties  as 
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from  which  we  get  the  total  lift  coefficient 
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Because  of  the  nature  of  the  pressure  loading  function,  equation  (54), 
equations  (66),  (67)  and  (68)  must  be  integrated  numerically. 

A  similar  analysis  leads  to  the  moment  equation  about  the  y  axis 
(see  Figs.  1  and  2)  as 
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'  /4P„  *  d4„  +  /  APE  X  dAE 
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For  the  first  integral  in  eq.  (69),  we  note  that 
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*LE 


+  Cc 
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and 

dAw  =  dxdy  =  (CwdC)(|  dn)  (71) 

and 

=  y  tan  A  =  n  tanA 


(72) 
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Using  these  definitions,  the  first  of  these  integrals  over  the  wing  can 
be  written  as 


0-1 


ll  =  qb  f  {|otanALCw  f  ACp  dig  +  Cw2  f  ACp  €wd?}dn  (73) 
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In  a  similar  manner,  the  second  of  these  integrals  become 
0-1  1  1 
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The  moment  coefficient  about  t.he  y  axis  is  defined  as 


cM  =  m  / (qSc) 
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In  order  to  evaluate  eqs.  (66),  (67),  (68),  and  (75),  we  must  first 
determine  the  constants,  in  eq.  (54)  for  both  the  wing  and  elevon. 

Since  the  elevon  solution  cannot  affect  the  wing  solution  (zone  of  silence). 


then  the  wing  loading  is  solved  independently  of  any  elevon  considerations 
through  a  Gaussian  reduction  of  eq.  (59). 

For  the  elevon  solution,  we  must  consider  the  downwash  produced  by 
that  portion  of  the  wing  in  the  forecone  of  some  arbitrary  point  on  the 
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elevon.  This  downwash  from  the  wing  on  the  elevon  will  enter  the  elevon 
solution  through  an  additional  term  in  the  boundary  conditions.  That 
is,  for  the  elevon 
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i=l 


(AK^tACp  «50.y0>]± 
E 
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where  w  is  the  downwash  produced  by  the  wing  at  a  particular  control 
w 

point  on  the  elevon.  Evaluation  of  w  is  considered  in  the  next  section. 

w 

Wing  Downwash  on  Elevon 

To  evaluate  the  downwash  from  the  wing  on  the  elevon,  a  direct 
approach  was  first  considered.  That  is,  integration  was  to  be  carried 
out  over  the  forecone  from  a  point  on  the  elevon  surface.  However,  because 
of  the  gap  between  wing  and  elevon,  this  did  not  seem  to  be  a  feasible 
approach,  especially  in  light  of  its  complexity.  Thus  another  approach  was 
considered  and  used. 

Since  the  wing  solution  is  completed  first,  sectional  properties  on 
the  wing  may  be  computed  at  any  spanwise  location.  Thus  the  wing  is 
divided  into  equally  spaced  "strips"  along  the  span  and  the  sectional  lift 
is  computed  for  each  strip.  It  is  assumed  that  a  horseshoe  vortex  is 
attached  to  each  strip  in  a  lattice  manner  whose  strength  is  determined  by 
the  sectional  lift  coefficient.  That  is,  from  the  Kutta- Joukowski  theorem, 
the  lift  per  unit  span  for  some  ith  section  is 


(AL/An)1  =  pvr  =  2qTi 


(77) 


vr 


i 


where 


Placement  of  these  vortices  are  shown  in  Fig.  20.  Note  that  in  the 


figure,  the  bound  portions  are  shown  attached  to  the  quarter  chord  and 
the  trailing  legs  trail  b< hind  the  wing  to  infinity. 

If  the  wing  is  symmetrically  loaded,  the  lift  may  be  written  as 
1 

L  =  qb  J  (cc£)w  dnw  (78) 
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so  that  for  a  strip 
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Dividing  by  An  and  equating  to  eq.  (77)  yields 


ri=‘^  / 


(cc.)  dn  ] . 
£  w  w  i 


If  it  is  assumed  that  over  th.  strip  the  term  cc;  is  constant,  then  the 
integration  may  be  completed  in  closed  form. 

At  the  junction  between  any  two  strips,  the  resulting  trailing  fila¬ 
ment  strength  is  the  difference  between  adjacent  filaments  and  may  be 
written  as 


(cC£)i+l 


(79) 


The  starting  point  of  these  filaments  is  assumed  to  be  the  wing 
trailing  edge  (for  lack  of  a  better  assumption).  Fortunately,  the  down- 
wash  induced  from  all  trailing  filaments  is  small  in  comparison  to  the 


I 
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free  stream  components  and  thus  an  error  in  the  starting  point  of  the 
filament  does  not  significantly  influence  the  final  elevon  loading. 
Obviously,  though,  the  filaments  must  start  somewhere  on  the  wing.  From 
From  a  planform  view  of  the  wing-elevon  trailing  filament  configuration, 
it  may  be  shown  that  if  a  particular  field  point  (x,y,z)  is  in  the  region 
of  influence  of  a  particular  filament  starting  at  (xp,yp),  then 

x  >  (xp+8 | y^-r | )  (80) 

where  r2  =  y2  +  z2.  That  is,  if  eq.  (80)  is  satisfied,  then  the  field 
point  (x,y,z)  "feels"  the  downwash  from  the  trailing  filament  (xp.y^). 

To  compute  the  velocity  induced  at  some  field  point  (x,y,z),  the 
Biot-Savart  Law  is  used  as  an  approximation.  Although  this  theorem  is 
actually  only  valid  for  subsonic  incompressible  flow,  it  will  be  suf¬ 
ficient  for  the  present  analysis  if  proper  account  is  taken  for  the 
region  of  influence  (zone  of  silence). 

In  vector  form,  the  Biot-Savart  relation  governing  the  velocity 
induced  at  some  point  by  a  vortex  filament  is 


V(x,y,z)  =  ~  ^  f 


r  x  di 


In  this  expression,  r  is  the  vector  from  the  field  point  to  an 
arbitrary  point  on  the  filament,  d?,  is  a  vectored  differential  length 
along  the  filament,  r  is  the  magnitude  of  the  w.;cor  r,  and  F  is  the 
filament  strength.  Referring  to  Fig.  21,  thpse  quantities  may  be 
written  as 


d£  =  (i  cos0p  +  k  sin0p)d£ 
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and 

T  =  (xr-x+L  cos6r)i  +  (y;-y)j  +  <Vln0r  '  z)k  (83) 

Equations  (82)  and  (83)  are  introduced  into  eq.  (81)  so  that  the 
integrals  for  the  three  velocity  components  are 
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These  integrals  may  be  evaluated  in  closed  form  with  the  results 


=  -T  (y^sinQ^,)  1/  (4tt) 

(92) 

=  f (x^sinGp  +  z  cosOp) 1/ (4n) 

(93) 

and 

Vz  =  T  (y^cosGp)  1/  (4tt)  (94 ) 

The  induced  velocity  from  these  filaments  on  the  elevon  is  of 
primary  concern  and  more  specifically  the  velocity  induced  normal  to  the 
elevon.  In  vector  form,  a  unit  normal  to  the  elevon  is 

n  =  -i  sin6+k  cosS  (95) 

so  that 

V  =  V*n  =  -  V  simS-V  cosd  (96) 

n  x  z 

or 

V  =  -Fy  I[cos(6r+6)  ]  /  (4tt)  (97  ) 

n  I  i 

Wing  Trailing  Edge  Blocking  of  Elevon 

To  complete  the  theoretical  analysis,  the  wing  thickness  and  blunt 
wing  trailing  edge  and  the  subsequent  effect  on  the  elevon  loading  must 
be  considered  in  addition  to  the  gap  between  the  wing  trailing  edge  and 
the  elevon  leading  edge.  The  con i igurat ion  being  used  in  this  analysis 
is  shown  in  Figure  22-  The  flow  just  aft  ot  the  blunt  wing  trailing 
edge  is  highly  turbulent  as  expected  and  details  describing  the  sub¬ 
sequent  flow  field  is  fundamentally  grounded  in  solutions  of  the  Navier 
Stokes  equations.  The  implication  then  is  that  viscous  effects  doming 
the  solution  as  was  suggested  in  Ref.^  However,  this  conclusion  is 

not  necessarily  true.  It  is  clearly  recognized  that  viscous  effects 


immediately  behind  the  wing  trailing  edge  are  very  important,  but  it  is 
believed  that  the  dominant  effects  over  most  of  the  elevon  can  be  pre¬ 
dicted  by  potential  flow  methods.  Justification  of  this  conclusion  is 
presented  in  a  later  section. 

As  depicted  in  Fig.  22,  the  flow  is  assumed  to  be  tangent  to  the 
surface  of  the  aft  portion  of  the  wing  for  all  angles  of  attack.  The 
streamlines  thus  emanating  from  the  upper  and  lower  aft  wing  surface 
impinge  on  the  elevon  surface  as  shown  in  Jig.  22,  forming  a  "blanked  out 
region  bounded  by  the  blunt  wing  trailing  edge,  the  upper  and  lower 
streamlines  and  a  section  of  the  elevon  chord  line.  In  this  blanked 
region,  details  of  the  flow  field  are  not  well  known.  Thus  for  purposes 
of  this  analysis,  it  is  assumed  that  this  wake  pressure  coefficient  is 
zero  such  that  essentially  free  stream  conditions  exist. 

From  Fig.  22  the  location  of  the  aft  impingement  point  is 
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That  is,  for  the  schematic  in  K  g.  2i,  ^  is  about 


(99) 

portion  of  the  wing. 
2.0"  and  1  is  about 


-10'. 


Equations  (98)  and  (99)  are  valid  tor  ‘  0.0.  If  the  elevon  defle 

tion  angle  is  positive,  then  the  forward  location  is  determined  from 
eq.  (98)  and  the  aft  location  from  eq.  (99).  It  is  obvious  from  Fig.  22 
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that  there  are  elevon  angles  where  the  leading  edge  of  the  elevon  lies 
entirely  within  the  wake  region  with  some  portion  of  the  aft  section  out¬ 
side  the  wake.  Also,  for  very  small  elevon  deflection  angles,  the  entire 
elevon  is  immersed  in  the  wing  wake  and  thus  is  completely  unloaded. 

In  order  to  justify  the  use  of  eqs.  (98)  and  (99)  to  form  a  blanked- 
out  region  (Newtonian  approach)  a  wind  tunnel  experiment  was  designed  to 
study  this  flow  field,  a  two-dimensional  model  consisting  of  a  wing  with 
a  blunt  trailing  edge  and  a  diamond  shaped  elevon  was  constructed  (see 
Fig.  22)  and  tested  at  a  Mach  number  of  1.94  in  a  4"  by  4"  test  section. 
The  total  pressure  was  about  30  psia  with  a  freest  ream  unit  Reynolds 
number  of  1.6  x  106.  Selected  schlieren  photographs  of  the  results  of 
these  tests  are  shown  in  Fig.  23.  In  the  upper  left  where  a  =  0.0°  and 
6  =  0.0°,  the  elevon  lies  entirelv  within  the  wake  of  the  wing  as  is 
clearly  seen.  The  important  thing  to  note  here,  though,  is  the  well  de¬ 
fined  streamline  emanating  I rom  the  top  surface  of  the  wing  and  impinging 
on  the  elevon  boundary  layer  near  the  apex  of  the  elevon  midchord.  In 
the  upper  right  where  i  =  0.(1'  and  ‘  -  -‘*.0",  the  leading  e.ge  of  the 
elevon  is  about  4°  outside  tin  wing  wake  region  hut  apparently  still 
inside  the  wing  boundary  1  aver.  I  lie  at  t  impingement  point  on  the  upper 
elevon  surface  is  still  ileutlv  isible  hut  the  l^rward  point  i  aiinut  tie 
clearly  located  . 

In  the  lower  left,  where  ■  =  -4.14  ami  =  --.O' ,  the  aft  impingement 
point  is  located  easily  while  the  elevon  leading  edge  remains  submerged  In 
the  wing  wake.  In  the  lower  right  where  <  =  -4.14”  and  =  -10.  (I0,  the 
results  are  inconclusive.  The  overall  results  of  these  tests  indicate 
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that  some  sort  of  blanking  region  is  formed  which  is  not  totally  inde¬ 
pendent  of  angle  of  attack  and  is  more  strongly  dependent  on  the  elevon 
deflection  angle  than  that  predicted  by  eqs.  (98)  and  (99).  Nevertheless 
these  equations  do  provide  "bounds"  for  the  mathematical  analysis  and,  as 
will  be  shown  later,  provide  good  agreement  with  experimental  data. 


VII.  SUPERSONIC  FLOW  -  RESULTS 


Comparison  with  experimental  data  of  the  theoretical  analysis  as 
previously  presented  is  divided  into  two  major  categories:  (1)  planar 
wings  and  (2)  deformed  wings.  In  many  cases,  the  planform  shapes  chosen 
were  dictated  by  available  experimental  data.  In  some  cases,  experimental 
data  was  available  for  both  deformed  and  undeformed  wings  with  the  same 
planform  shape.  Four  of  the  planform  shapes  analyzed  are  shown  in  Fig.  24. 

Planar  Wings 

Results  of  the  pressure  loading  and  sectional  lift  and  pitching 

moment  for  three  of  these  configurations  are  shown  in  Figs.  25-30  for 

selected  Mach  numbers  and  low  angles  of  attack.  More  complete  data  summarv 

results  for  these  configurations  are  presented  in  Appendix  A.  tn  these 

figures  the  theoretical  solution  for  the  pressure  loading  is  basic  three- 

dimensional  supersonic  theory  as  presented  in  Refs.  13  and  18.  No 

vorticity  paneling  is  required.  The  double  summ.it  ton  terms  in  eq .  (54)  are 

omitted  and  only  the  P(x  ,y  )  sin  (i)  term  is  used  to  compute  the  pressure. 

o  o 

In  most  cases,  the  sectional  normal  forces  agree  well  with  the  exper¬ 
imental  data  especially  on  the  inboard  sections  ot  the  wings.  It  is 
observed  that  the  worst  agreement  is  near  the  wing  tips  as  one  would 
expect.  This  deviation  at  the  tips  is  not  as  detrimental  as  might  first 
appear  since  most  of  the  configuration  loading  is  concentrated  on  the 
inboard  sections. 
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The  sectional  moment  coefficients  about  the  y-axis  shows  the  same 
kind  of  agreement  as  does  the  normal  forces.  In  each  case  the  dis¬ 
crepancies  are  not  severe  and  the  location  of  the  sectional  aerodynamic 
center  is  predicted  quite  well. 

The  experimental  data  for  these  plots  (including  those  in  Appendix 
A)  were  obtained  from  Refs.  29  and  30  and  was  integrated  numerically  to 
obtain  the  sectional  normal  forces  and  pitching  moments.  Integration  to 
produce  the  normal  forces  agreed  quite  well  with  that  tabulated  in  the 
original  reports;  however,  the  sectional  pitching  moment  data  as  pub¬ 
lished  in  Refs.  29  and  30  was  not  reproducible  by  numerical  means.  Con¬ 
sequently,  these  experimental  data  as  presented  in  Figs.  25-30  (and 
Appendix  A)  for  the  sectional  coefficients  was  obtained  by  numerical 
integration  of  the  pressures  and  is  not  that  which  is  "tabulated"  in 
Refs.  29  and  30. 

For  these  planar  configurations,  the  theory  is  a  linear  theory  and 
thus  tabulation  of  the  total  forces  and  moments  is  sufficient  to  illus¬ 
trate  overall  agreement.  Table  2  is  a  comparison  of  the  lift  curve  slope 
and  moment  curve  slope  for  the  three  planar  configurations  under  study, 
as  found  by  numerical  integration  of  the  curves  in  Figs.  26,  28  and  30. 

From  Figs.  25-30  and  those  in  Appendix  A,  it  is  difficult  to  picture 
exactly  how  the  sectional  loading  profiles  are  spread  over  the  wing. 
Consequently,  Figs.  31-34  were  produced  showing  the  isobars  and  the  large 
influence  produced  by  the  wing  tips.  Note  also  the  supersonic  leading  edge 
behavior  for  the  low  aspect  ratio  trapezoidal  wing  and  rectangular  wing. 

It  is  apparent  from  Figs.  25-30  and  from  Table  2  that  the  linear 
theory  is  excellent  for  most  cases  considered.  The  disagreement  near  the 
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Table  2 


Low  AR 

High  AR 

Delta 

Trapezoidal 

Trapezoidal 

Wing 

Wing 

Wing 

Experiment 

.029 

.033 

.048 

Theory 

.029 

.032 

.053 

,  Experiment  -.029  -.019  -.064 

M  'v 

a  y - — — - 

a  =  0  Theory  -.029  -.019  -.073 


leading  edge  and  wing  tip  is  to  be  expected  and  is  due  to  wing  thickness, 
shocks,  and  wing  twist;  the  latter,  of  course,  is  not  significant  for  the 
low  aspect  ratio  wings. 

Deformed  Wings 

These  same  wings  under  various  kinds  of  deformation  were  also 

analyzed  by  including  vorticity  paneling  as  well  as  the  3-D  theoretical 

loading  terms  in  the  AC  distribution.  Deformation  shapes  of  three  of 

P 

these  wings  are  shown  in  Fig.  35.  The  delta  wing  (No.  1)  in  Fig.  35 
was  cambered  and  the  slope  relative  to  the  chord  line  of  the  mean 
camber  line  was  determined  to  be 

0  =  -3.00825  +  tan-1 [tan  6°(f-£n+n)]  (100) 

where  n=y/(b/2)  and  £=(x-x  )/c.  The  camber  (see  Fig.  35)  is  uniform  in 

the  spanwise  direction  and  produces  a  maximum  deflection  angle  of  three 
degrees  relative  to  the  chord  line.  Although  this  seems  like  a  small 
task,  this  is  actually  a  severe  test  of  the  theory  and  its  agreement  with 
experimental  data.  In  Figs.  36  through  41,  resul ts  are  presented  for  three 
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Fig.  35.  Schematic  of  Three  Deformed  Wing  Shapes. 
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Figure  36.  Pressure  Distribution  on  a  Cambered  Delta  Wing;  M  =  2.01,  a 
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Figure  36.  Continued 


CCMY  CCN 

.140  -0.80  -1.20  -1.60  _0. 00  0.20  0.40 


94 


o  CflHBEREO  DELTA  WING 

ao 


o 


o 

CO 

□ 


ALPHA  =  6.0000 

LAMDA  =  70.0000 

MACH  NO.  =  2.0100 

CR  =  2. 7475 
CT  =  0.0000 


LOO  o'.  25  oVsO  ~CL75  l\  00 

Y/  (B/2) 


CAMBERED  0ELTR  WING 


ALPHA  =  6.0000 

LAMOA  =  70.0000 

MACH  NO.  =  2.0100 

CR  =  2.7475 
CT  =  0.0000 


o 

o 

^'.00  o'.  25  oT50  CL  7  b  l'.OO 

Y/  (B/2) 
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Figure  38.  Pressure  Distribution  on  a  Cambered  and  Twisted  Delta  Wing; 
M  =  2.01,  a  =  6.0. 
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Figure  41.  Sectional  Normal  Force  and  Pitching  Moment  for  a 
Warped  Trapezoidal  Wing;  M  =  2.01;  ot  =  6.0. 
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deformed  configurations  at  a  Mach  number  of  2.01  and  angle  of  attack  of 
6  degrees.  A  more  complete  summary  of  agreement  between  the  theory  and 
experimental  data  is  presented  in  Appendix  A. 

For  the  cambered  delta  wing,  agreement  at  the  inboard  stations  is 
relatively  good  while  outboard  agreement  is  rather  poor.  However,  note 
that  in  the  regions  where  most  of  the  wing  loading  occurs,  the  agreement 
is  satisfactory.  Additionally,  the  inclusion  of  the  doublet  paneling 
terms  in  the  mathematical  formulation  provides  the  correct  trends  and  an 
added  degree  of  accuracy. 

For  the  same  planform  shape,  linear  spanwise  twist  is  added  to  the 
camber  deformation  such  that  the  mean  camber  surface  is  now  defined  as 

0  .  =  tan-1 [0.140541(£-£n+n)n]  (101) 

1  >  J 

Results  of  this  configuration  are  shown  in  Figs.  38  and  39  and  in 
Appendix  A.  Similar  results  are  obtained  for  these  cases  as  was  observed 
for  the  cambered  (alone)  case.  Agreement  is  acceptable  at  the  inboard 
stations  but  is  questionable  at  the  outboard  stations. 

In  observing  the  experimental  data  for  both  the  cambered  alone  and 
the  cambered  and  twisted  wings,  the  outboard  stations  (r,  .5)  appear  to 

have  a  supersonic  leading  edge.  Whereas  the  high  leading  edge  sweep 
angle  gives  rise  to  a  subsonic  leading  edge.  In  addition,  at  the  n  =  0.5 
station,  a  weak  shock  appears  to  be  located  near  £  =  .4  for  M  =  1.61  but 
is  not  present  for  M  =  2.01  (see  Appendix  A).  Since  the  theory  does  not 
account  for  the  shocks,  the  disagreement  in  this  region  is  not  surprising. 

The  next  wing  subjected  to  the  theoretical  analysis  was  a  "warped" 
trapezoidal  wing  (see  Fig.  35).  For  both  M  =  1.61  and  2.01,  the  leading 


102 

edge  is  supersonic  and  multiple  "fold  over"  regions  occur  because  of  the 
low  aspect  ratio  of  the  fact  that  the  wing  is  trapezoidally  shaped.  The 
equation  governing  the  mean  warped  camber  surface  is 

©i  =  -tan-1  [-0.033454tt  sin(-^p)sin(n£;)  ]  (102) 

Results  for  this  wing  are  presented  in  Figs.  40  and  41  and  in 
Appendix  A.  Considering  the  wing  deformation,  the  agreement  is  remark¬ 
ably  good  over  the  entire  wing  for  both  Mach  numbers.  Note  also  that  the 
addition  of  vorticity  paneling  to  the  AC^  distribution  predicts  the  right 
trends  and  correctly  accounts  for  local  wing  deformations  in  the  upstream 


running  Mach  cones. 
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Elevon  Loading 

Very  little  experimental  data  exists  for  loading  on  an  elevon 

immersed  in  the  wake  of  a  wing  in  supersonic  flow.  The  only  data  found 

was  Nielsen  (Ref.  1)  which  included  the  blunt  wing  trailing  edge  as 

previously  discussed.  Figure  42  is  a  comparison  of  the  normal  force 

coefficient,  C_.,  and  hinge  moment,  C  ,  versus  wing  angle  of  attack  for 
N  H 

an  elevon  deflection  angle  of  -10°.  In  this  figure,  the  present  theo¬ 
retical  results  are  compared  to  experimental  data  and  to  the  theoretical 
results  from  Ref.  1.  In  the  present  theory,  the  blanking  effect  from 
the  wing  is  included  which  seems  to  properly  account  for  the  normal  force 
coefficient  variation.  It  is  noted  also  that  compressibility  or  the 
Mach  number  effect,  is  also  properly  accounted  for,  thereby  providing 
excellent  agreement  (at  least  for  this  case)  between  experimental  data 
and  theory. 

The  hinge  moment  for  the  elevon  is  also  plotted  and  agreement  here 
is  not  good  except  at  very  small  angles  of  attack.  Reasons  for  this 
discrepancy  may  be  the  result  of  one  of  the  following  problems  in  the 
analysis  or  perhaps  combinations  thereof:  (1)  It  is  known  that  viscous 
flow  in  this  region  could  cause  large  changes  in  the  pressure  field  as 
originally  hypothesized  in  Ref.  1.  However,  because  of  the  excellent 
agreement  of  the  normal  force  coefficient  data  from  a  potential  flow 
analysis  and  the  schlieren  photographs  of  this  region,  it  is  certainly 
not  conclusive  that  viscous  effects  dominate  the  solution.  (2)  The 
downwash  model  as  outlined  in  the  theory  section  is  inadequate  in  pre¬ 
dicting  the  correct  downwash  velocities.  As  a  matter  of  fact,  if  the 
downwash  distribution  is  assumed  to  be  constant  across  the  elevon  span 
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and  chord  slightly  better  results  are  obtained.  However,  since  the  total 
normal  force  is  predicted  quite  well,  the  chordwise  distribution  is 
apparently  in  error.  From  the  schlieren  photographs  (Fig.  23)  it  appears 
that  much  more  of  the  upper  surface  of  the  elevon  may  be  blanked  out  than 
that  predicted  in  the  theory.  If  this  were  true,  the  normal  force  coef¬ 
ficients  would  decrease  slightly  and  the  moment  coefficients  would 
increase  (leading  edge  up)  which  would  give  rise  to  an  increase  in  the 
associated  hinge  moment.  (3)  The  assumed  location  of  the  trailing 
vortices  relative  to  the  elevon  is  perhaps  incorrect.  (4)  In  analyzing 
the  experimental  data  from  Ref.  1,  there  is  possibly  a  problem  in  inter¬ 
pretation  and  conversion  to  the  reference  areas  and  lengths  used  in  the 
present  theoretical  analysis. 

In  conclusion,  it  is  recognized  that  the  flow  field  immediately 
behind  the  blunt  wing  trailing  edge  and  elevon  leading  edge  is  not  well 
understood.  Detailed  experiments  to  measure  pressure  and  loads  accom¬ 
panied  by  schlieren  photographs  are  required  before  an  adequate  model  can 
be  postulated. 

Finally,  Figure  43  is  a  plot  of  the  theoretically  predicted  normal 
force  and  hinge  moments  for  the  elevon  when  a  gap  exists  between  the 
elevon  and  the  wing.  Similar  results  are  obtained  as  for  the  "no  gap" 
case  but  the  influence  of  the  gap  is  clearly  identified  and  the  present 
theory  predicts  the  correct  trends. 
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Figure  A7.  Pressure  Distribution  on  a  Flat  Trapezoidal  Wing;  M  =  2.01, 
a  =  4. 
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Figure  All.  Pressure  Distribution  on  a  Flat  Trapezoidal  Wing;  M  =  2.01, 

a  =  10.0. 
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Figure  A14.  Sectional  Normal  Force  and  Pitching  Moment  for  a 
Flat  Delta  Wing;  M  =  1.61,  a  =  4.0. 
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Pressure  Distribution  on  a  Flat  Delta  Wing;  M  =  1.61,  01  =  10.0. 
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Figure  A22.  Sectional  Normal  Force  and  Pitching  Moment  for  a 
Flat  Tapered  Wing;  M  =  1.61,  a  =  2.0. 
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Figure  A23.  Pressure  Distribution  on  a  Flat  Tapered  Wing,  M  =  2.01,  a  =  2.0 
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Figure  A25.  Continued 
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Figure  A26.  Sectional  Normal  Force  and  Pitching  Moment  for 
a  Flat  Tapered  Wing;  M  =  1.61,  a  =  4.0. 
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Figure  A27.  Continued 
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Figure  A28.  Sectional  Normal  Force  and  Pitching  Moment  for  a 
Flat  Tapered  Wing;  M  =  1.61,  a  =  10.0. 
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Figure  A29.  Continued 
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Figure  A30.  Sectional  Normal  Force  and  Pitching  Moment  for  a 
Cambered  Delta  Wing;  M  =  1.61,  a  =  6.0. 
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Figure  A31.  Pressure  Distribution  on  a  Cambered  Delta  Wing;  M  =  2.01,  a 
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Figure  A31.  Continued 
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Figure  A32.  Sectional  Normal  Force  and  Pitching  Moment  for  a 
Cambered  Delta  Wing;  M  =  2.01,  a  =  6.0. 
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Figure  A33.  Pressure  Distribution  on  a  Cambered  and  Twisted  Delta  Wing. 
M  =  1.61,  a  =  6.0. 
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Figure  A3A . 


Sectional  Normal  Force  and  Pitching  Moment  for  a 
Cambered  and  Twisted  Delta  Wing;  M  -  1.61,  a  =  6.0. 
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Figure  A35 . 


Pressure  Distribution  on  a  Cambered  and  Twisted  Delta  Wing; 
M  =  2.01,  a  =  6.0. 
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Figure  A35.  Continued 
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Figure  AJ6 .  Sectional  Normal  Force  and  Pitching  Moment  for  a 
Cambered  and  Twisted  Delta  Wing;  M  =  2.01,  a  =  i 
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Figure  A37.  Pressure  Distribution  on  a  Warped  Trapezoidal  Wing;  M  =  1.61, 

=  6.0. 
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Figure  A37.  Continued 
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Figure  A39 .  Pressure  Distribution  on  a  Warped  Trapezoidal  Wing;  M  =  2.01 
a  =  6.0. 
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Figure  A39 .  Continued 


SUBSONIC  COMPUTER  PROGRAM  USERS  GUIDE 

I.  General  Description  of  the  Main  Program  and  Subroutines  for 
the  Subsonic  Program 

MAIN  PROGRAM  -  This  is  the  executive  segment  which  primarily  controls  and 
coordinates  the  various  tasks  performed  by  the  subprograms.  Other 
functions  performed  by  the  main  program  include:  (a)  all  input  and  the 
majority  of  output;  (b)  non-dimensionalization;  (c)  calculation  of  the 
boundary  conditions  (the  "B"  matrix);  (d)  computation  of  all  total 
characteristics,  both  exact  results  and  numerical  integrations  (this 
includes  wing-elevon,  fuselage,  suction,  and  all  combinations  of  these). 

BLOCK  DATA  -  This  data  gives  coordinates  of  a  non-planar  mean  line  (XTB.ZTB) 
as  well  as  the  slope  (ZPTB)  for  use  in  calculating  wing  deformations. 

The  present  data  is  for  the  deformed  Hawk  wing. 

SUBROUTINE  SURFIN  -  This  routine  computes  the  value  of  the  surface  integral 
(the  "A"  matrix)  by  summing  the  effects  of  each  grid  panel  at  a  control 
point.  It  is  called  once  for  the  wing  and  once  for  the  elevon  for  each 
control  point. 

SUBROUTINE  SECLM  -  The  sectional  loading  characteristics  are  calculated  for 
the  wing  or  elevon  in  this  subroutine.  It  is  called  after  the  pressure 
loading  coefficients  are  determined. 

SUBROUTINE  PRESUR  -  This  subroutine  is  used,  after  the  pressure  loading 
coefficients  are  determined,  to  calculate  and  print  AC^  at  NSPAN  by 
NCORD  points  on  the  lifting  surfaces. 
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SUBROUTINE  GEOM  -  This  subroutine  sets  up  the  geometry  of  both  the  grid 
panels  and  the  control  points,  and  then  centers  the  control  points 
inside  the  grid  element  in  which  they  appear.  Basically,  a  cosine 
distribution  is  used  for  both  the  grid  elements  and  control  points. 

This  routine  is  called  once  for  the  wing  and  once  for  the  elevon. 

SUBROUTINE  INTV  -  This  routine  contains  the  actual  equations  of  the  sur¬ 
face  integral  and,  therefore,  is  utilized  by  SUBROUTINE  SURFIN. 

SUBROUTINE  INVERT  -  Gauss  elimination  is  used  in  this  subroutine  to 

"invert"  the  A  matrix  and  solve  for  the  unknown  leading  coefficients, 

B  's. 
nm 

SUBROUTINE  HINGE  -  This  contains  the  equations  of  the  closed  form  solution 
for  the  elevon  hinge  moment. 

SUBROUTINE  TLUP  -  This  is  a  "table-lookup"  routine  which  uses  the  data 
given  by  the  variables  in  BLOCK  DATA  to  interpolate  for  wing  mean  line 
slopes  used  in  determining  deformation  angles.  It  is  called  by 
FUNCTION  WARP. 

FUNCTION  WARP  -  This  function  uses  either  the  table-lookup  or  a  function 
definition  to  find  wing/elevon  deformation.  It  is  called  by  both  the 
main  program  to  determine  deformation  of  control  points  (included  in 
the  boundary  conditions)  and  by  SUBROUTINE  SURFIN  to  determine  deforma¬ 
tion  of  grid  panels. 
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II.  INPUT  DATA 


A.  Data  Cards: 

Card  No. 

Format 

Variable  List 

1 

18A4 

WING 

2 

7F10.0 

SP AN , CR , CT , LAMDA , S AREA , CBAR , MACH 

3 

FIO.0,15 

ALPHA, KSTOP 

A 

1015 

NC , NS , NX , NY , NCF , NSF , NXF , NYF , NSPAN 
NCORD 

5 

5F10.0 

CRF , CTF , DELTA , ALTWMX , EPS 

6 

F10 . 0, 15 , 2F10 . 0 

HNL , IDEF , AMP , DB 

(Additional 

cards) 

3F10.0 

ALPHA , DELTA , ALTWMX 

("Additional  cards"  allows  more  than  one  case  of  the  given 
geometry  to  be  run  with  different  boundary  conditions. 
There  should  be  (KSTOP-1)  additional  cards.) 


B.  Variable  definitions  (all  linear  dimension  units  must  be  consistent) 
WING  -  Alphanumeric  problem  description 
SPAN  -  Total  wing  span 
CR  -  Wing  root  chord 
CT  -  Wing  tip  chord 

LAMDA  -  Wing  leading  edge  sweep  angle 

SAREA,CBAR  -  Reference  area  and  reference  length.  If  0.0  is  input, 
the  total  area  and  mean  aerodynamic  chord  will  be 
calculated  and  used. 

MACH  -  Mach  number 

ALPHA  -  Angles  of  attack  of  wing,  in  degrees.  Leading  edge  up 
is  positive. 

KSTOP  -  Number  of  cases  of  different  boundary  conditions  that  are 
to  be  run. 
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NC,NS  -  Number  of  chordwise  and  spanwise  control  points  (re¬ 
spectively)  on  the  wing  (typically  4x4) . 

NX, NY  -  Number  of  chordwise  and  spanwise  (respectively)  grid  panels 
on  the  wing  (typically  45x45). 

NCF.NSF  -  Chordwise  and  spanwise  control  points  on  the  elevon 
(typically  3x4).  NSF  must  be  equal  to  NS. 

NXF,NYF  -  Chordwise  and  spanwise  grid  panels  on  the  elevon 
(typically  20x45).  NYF  must  be  equal  to  NY. 

NSPAN,NCORD  -  Number  of  spanwise  by  chordwise  points  at  which  AC 
distribution  is  to  be  printed  out.  P 

CRF  -  Elevon  root  chord 

CTF  -  Elevon  tip  chord 

DELTA  -  Elevon  deflection  angle,  in  degrees.  Trailing  edge  downward 
is  positive. 

ALTWAX  -  Maximum  twist  angle  of  the  wing  tip,  in  degrees.  Leading 
edge  down  at  the  tip  is  positive. 

EPS  -  Gap  distance  between  wing  and  elevon. 

HNI  -  Hinge  line  location  from  elevon  leading  edge. 

IDEF  -  Deformation  check.  If  IDEF=1,  the  Table-Lookup  Routine  is 
used.  If  IDEF=0,  a  function  definition  is  used. 

AMP  -  Amplitude  of  a  function  deformation. 

DB  -  Diameter  of  fuselage. 


III.  OUTPUT  DATA 


Major  Headings  Are: 

A.  INPUT  DATA  -  Dimensional  imput  as  it  was  given  on  the  data  card 

B.  ACTUAL  VALUES  USED  IN  PROGRAM  -  Variables  non-dimensionalizci 

by  b/2 . 

C.  CONTROL  POINTS  -  Location  based  on  non-dlmensionalizat ion  of 
local  wing  chord  by  b/2.  Gap  distance  is  not  included  in  XO. 


AD-A101  675 


UNCLASSIFIED 

3"  3 
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D.  BOUNDARY  CONDITIONS 

I  -  Control  point  number. 

B  -  Boundary  condition  matrix  values. 

XO,YO  -  Control  point  location  as  above,  except  with  gap  distance 
included. 

ALTW  -  Local  twist  angle  at  the  control  point  (degrees). 

TAU  -  Local  deformation  angle  at  the  control  point  (degrees). 

E.  POLYNOMIAL  COEFFICIENTS  -  B  ' s  of  pressure  loading  function. 

nm 

F.  PRESSURE  DISTRIBUTION 

ELEVON  LEADING  EDGE  -  Location  in  percent  local  chord  (gap  included) 

Y  -  Nondimensional  spanwise  location  of  the  given  pressure. 

X  -  Nondimensional  chordwise  location  in  percent  local  chord 
(gap  included) . 

XDIS  -  Chordwise  location  in  terms  of  nondimensionalization  by  b/2. 
DCP  -  ACp  (pressure  loading) 

C  -  Local  chord  length  of  the  wing  only. 

G.  SECTIONAL  CHARACTERISTICS 

CLW  -  cc^  of  wing  only. 

CMW  -  czc  of  wing  only 
m 

CLE  -  cc^  of  elevon  only 

CMH  -  c2c  of  elevon  only 
m 

CCL  -  cc^  of  total  wing-elevon. 

CCM  -  c2c  of  total  wing-elevon 
m 

XCP  -  Center  of  pressure  in  percent  local  chord  behind  the  leading 
edge. 

Y  -  Nondimensional  spanwise  location  of  the  sectional  characteristics. 
C  -  Local  wing  chord  nondimensionalized  by  b/2. 
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CF  -  Local  elevon  chord  nondimensionalized  by  b/2. 

CTDT  -  Total  local  wing-elevon  chord  (gap  not  included) 

CDW  -  cc_  of  wing  only 
1 

CDE  -  ccQ  of  elevon  only 

CCD  -  ccQ  of  total  wing-elevon 

H.  TOTAL  FORCES  AND  MOMENTS 

TOTAL  LIFT,  etc-  -  Characteristics  of  wing-elevon  combination. 

CORRECTED  LIFT  -  Lift  multiplied  by  sin  a. 

LIFT  AND  MOMENT  DUE  TO  SUCTION  -  Contribution  of  suction  only. 

Lift  value  is  correct  only  when  a>0.  Both  terms  are  incorrect 
for  6^0. 

FUSELAGE  LIFT  AND  MOMENT  -  Contribution  of  fuselage  only. 

WING- FUSELAGE. . .  -  Wing-elevon  fuselage  characteristics. 

TOTAL  CONFIGURATION...-  Wing-elevon-f uselage-suction  characteristics. 

HINGE  MOMENT  ABOUT  ELEVON  HINGE  LINE  -  Elevon  hinge  moment. 

IV.  NON-PLANAR  CONFIGURATION  OPTIONS 

A.  Wing  Twist  -  Statement  numbers  2310  and  4540  calculate  wing  twist 

with  a  linear  distribution  along  the  span  (Y0(J)).  Other  (spanwise) 

variations  may  be  used  by  substituting  the  proper  equation. 

B.  Deformations 

1.  Statement  numbers  2200  and  5430  must  be  deleted  if  the  elevon 
is  to  be  deformed  in  addition  to  the  wing. 

2.  The  FUNCTION  WARP  routine  supplies  the  deformation  angles 

a.  If  the  input  variable  IDEF  is  0,  the  equation  for  WARP 
given  by  statement  number  9430  is  used  to  compute  deforma¬ 
tions.  The  variable  (also  input)  AMP  may  be  used  to 
modify  this  equation.  The  supplied  equation  must  be  for 
the  slope  (dz/d£)  of  the  mean  line.  Different  functions 
could  be  used  at  different  spanwise  locations  if  the 
variables  YO(J)  in  statement  number  2160  and  YA  in  state¬ 
ment  number  5450  are  included  in  the  argument  list,  and 
the  appropriate  modifications  are  made  to  FUNCTION  WARP. 

b.  If  IDEF  >  1,  a  table  lookup  routine  is  used  to  compute  de¬ 
formations  based  on  a  numerical  table  given  by  the  variables 
in  BLOCK  DATA. 


SUPERSONIC  COMPUTER  PROGRAM  USERS  ''UIDE 

I.  General  Description  of  the  Main  Program  and  Subroutines 
for  the  Supersonic  Program 

MAIN  PROGRAM 

This  is  the  executive  portion  of  the  program  and  basically  controls 
the  sequencing  of  the  other  subroutines.  All  data  to  be  read  into  memory 
is  completed  here  as  well  as  all  nondimensionalization.  Decisions  are 
made  on  the  type  of  pressure  solutions  desired  and  appropriate  transfer 
of  control  is  made.  Geometry  for  the  vorticity  paneling  for  the  wing 
and  elevon  is  completed.  Wing  pressure  solutions  are  called  and  the 
trailing  legs  vorticity  strengths  are  computed.  Most  of  the  entire 
elevon  solution  is  complete  in  this  program  segment  and  appropriate 
computational  output  is  generated. 

Subroutine  BLOCK  DATA 

This  subroutine  establishes  the  spanwise  locations  at  which  the 
chordwise  pressure  distribution  si  computed.  If  one  wishes  to  alter  these 
locations,  it  must  be  done  in  this  subroutine. 

Subroutine  DAUGHN 

This  subroutine  calculates  the  downwash  of  any  arbitrary  location  on 
the  elevon  due  to  the  trailing  vortices  on  the  wing.  This  routine  is 
called  after  the  wing  loading  and  trailing  vorticity  strength  is  estab¬ 
lished. 
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Subroutine  EXACT 

This  subroutine  calculates  the  pressure  distributions  on  the  wing  or 
elevon  using  planar  three-dimensional  supersonic  flow  theory.  This 
subroutine  calls  the  load  coefficient  subroutine  where  sectional  forces 
and  moments  are  calculated. 

Subroutine  LODCOF 

This  subroutine  takes  the  AC^  arrays  and  integrates  the  chordwise 
distribution  using  Simpson's  Rule.  It  integrates  the  spanwist  distribu¬ 
tions  using  trapezoidal  integrations. 

Subroutine  DOWN 

This  subroutine  calculates  a  numerical  value  for  the  downwash  on 
the  wing  only  at  some  (x,y)  location.  It  is  used  primarily  to  check 
boundary  condition  agreement. 

Subroutine  SURF 

This  subroutine  computes  the  surface  integral  over  the  forecone 
from  the  field  point  (x,y)  by  summing  the  contribution  from  each  subpanel. 
In  doing  so  it  loads  the  "A"  matrix  to  be  inverted  later  for  the  unknown 
polynomial  coefficients,  B^. 

Function  CW 

This  function  evaluates  a  cosine  function  which  is  used  frequently 
in  the  program. 

Subroutine  INTV 

This  subroutine  evaluates  the  results  of  the  downwash  integral  over 


a  subpanel. 


I 
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Subroutine  INVERT 

This  subroutine  is  a  matrix  "inversion"  routine  which  actually  uses 
a  Gaussian  reduction  to  find  the  unknown  coefficients,  Bj  • 

Function  PRESS 

This  function  evaluates  the  appropriate  expression  to  find  the  AC^ 
at  a  point  using  supersonic  three-dimensional  theory. 

Subroutine  ELLIP 

This  subroutine  numerically  evaluates  an  elliptic  integral  used  in 
Function  Press. 

Subroutine  THA 

This  subroutine  establishes  the  deformation  equations  of  the  deformed 
wings  and  returns  the  local  wing  slope  relative  to  the  wing  chord  line. 

Subroutine  SUM 

This  subroutine  performs  a  summation  step  used  often  in  Subroutine 

SURF. 


Subroutine  ETHA 


This  subroutine  is  identical  to  subroutine  THA  except  the  deformation 


equations  define  local  slopes  on  the  elevon  surface. 

II.  Input  Data 

A.  Data  Cards: 


Card  No. 

Format 

Variable  List 

1 

815 

IDOWN,  IWGEX, IWGVP, IELEX, IELVP 
IPRINT,  IWGDF,  IELDF 

2 

6I5.2F10.0 

NC , NS , NX , NY , NCORD , NSPAN , VAG , DB 

3 

12A4 

CONFIG 

4 

4F10.0 

SPAN,CR,CT,LAMDA 

5 

4F10.0 

SWING, CBRW, MACH , ALPHA 

6 

415 , F10 . 0 

LVON ,NCE , NXE , NCRE , CRE 

7 

5F10.0 

CTE , GAP , DELTA , XH , S I V 

8 

2F10.0 

THIK,THTED 

1 
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Additional  configurations  can  be  run  by  adding  sets  of  cards 
4  through  8.  After  the  last  set  to  be  run,  SPAN  should  be 
read  in  as  0.0. 

B.  Variable  Definition  (all  linear  dimension  units  must  be  consistent) 

IF  (IDOWN  =  1)  a  subroutine  is  called  in  which  the  downwash  is 
computed  at  all  pressure  computed  points  on  the  wing. 

IF  (IWGEX  =  1)  uses  exact  solution  for  wing. 

IF  (IWGVP  =  1)  uses  vortex  paneling  also  for  wing  solution. 

IF  (IELEX  =  1)  uses  exact  solution  for  elevon. 

IF  (IELVP  =  1)  uses  vortex  paneling  also  for  solution  of  elevon. 

IF  (IPRINT  =  1)  a  more  extensive  printout  is  produced  at  inter¬ 
mediate  computation  points. 

IF  (IWGDF  =  1)  the  wing  has  deformation  as  specified  in  function 
THA.  Otherwise  wing  deformation  is  zero  (planar) . 

IF  (IELDF  =  1)  the  elevon  has  deformation  as  specified  in  function 
ETHA.  Otherwise  elevon  deformation  is  zero  (planar). 

NC  is  the  number  of  chordwise  control  points  for  the  wing. 

NS  is  the  number  of  spanwise  control  points  for  the  wing. 

NX  is  the  number  of  grid  points  in  the  chordwise  direction. 

NY  is  the  number  of  gridpoints  in  the  spanwise  direction. 

NCORD  is  the  number  of  chordwise  points  at  which  the  pressure  is 
computed  for  the  wing. 

NSPAN  is  the  number  of  spanwise  points  at  which  the  pressure  is 
computed  for  the  wing. 

VAG  is  the  angle  in  degrees  that  the  vortices  leave  the  wing 
reference  plane. 

DB  is  the  dimensional  body  diameter. 

CONFIG  is  the  alphanumeric  configuration  description. 

SPAN  is  the  dimensional  total  wing  span. 


CR  is  the  exposed  wing  root  chord. 

CT  is  the  wing  tip  chord. 

LAMDA  is  the  wing  leading  edge  sweep  angle  in  degrees. 

SWING  is  the  dimensional  reference  area.  If  0.0  is  read  in,  the 
program  will  compute  and  use  the  wing  area  as  the  reference 
area. 

CBRW  is  the  dimensional  reference  length.  If  0.0  is  read  in, 
the  program  will  compute  and  use  the  mean  aerodynamic  chord 
as  the  reference  length. 

MACH  is  the  freestream  Mach  number. 

ALPHA  is  the  angle  of  attack  of  the  wing  alone  in  degrees. 

IF(LVON  =  0)  the  elevon  is  omitted  from  all  computations  and 
the  wing  alone  case  is  analyzed. 

NCE  is  the  number  of  control  point  locations  on  the  elevon  in 
the  chordwise  direction. 

NXE  is  the  number  of  subelements  on  the  elevon  in  the  chordwise 
direction. 

NCRE  is  the  number  of  points  in  the  chordwise  direction  on  the 
elevon  at  which  pressure  is  to  be  computed. 

CRE  is  the  elevon  root  chord. 

CTE  is  the  elevon  tip  chord. 

GAP  is  the  gap  distance  between  the  wing  TE  and  the  elevon  LE. 

DELTA  is  the  elevon  deflection  angle  in  degrees  relative  to  wing 
trailing  edge  down  is  positive. 

XH  is  the  dimensional  distance  from  the  Y  axis  to  the  elevon 
hinge  line. 

SIV  is  the  nondimens ional  chordwise  location  at  which  the 
wing  trailing  vortices  are  started. 

THIK  is  the  dimensional  wing  TE  thickness. 

THTED  is  the  (positive)  wing  surface  slope  at  the  trailing  edge 
in  degrees. 
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Output  Data 
Major  Headings  Are: 

A.  Input  Data  -  Dimensional  or  nondimensional  data  as  it  was 
given  on  the  data  cards. 

B.  Actual  Values  used  in  Program  -  Variables  after  non- 
dimensionalization  by  b/2. 

C.  Pressure  Distribution  - 

D.  Sectional  Properties  as  Follows: 


YA 

«L*> 

CA 

CCL 

CCMY 

XAC 

XCPA 

(n) 

(c) 

(CCJl) 

(c2c  ) 
m 

y 

(L) 

ac 

(X  ) 
ac 

E.  Total  load  coefficients  from  "exact"  theory 

CN  -  Total  normal  force  coefficient  for  wing  only. 

CLT  -  Total  load  coefficient. 

CMTY  -  Total  moment  coefficient. 

YBAR  -  Spanwlse  location  for  mean  aerodynamic  chord. 

SIAC  -  Chordwise  location  for  aerodynamic  center. 

CMC 4  -  Moment  coefficient  about  the  quarter  chord. 

CD-CDMIN  -  Induced  drag. 

SAREA  -  Reference  area. 

CBAR  -  Reference  length. 

F.  [For  deformed  wings,  the  following  is  also  output] 

1.  (a)  Control  point  location;  XO,  YO,  local  chord;  CA  , 

The  x  location  of  the  last  row  of  control  points  is  re¬ 
located  at  trailing  edge  of  wing  and  new  locations  are 
printed  out. 

(b)  XO  and  YO  are  the  nondimensional  control  point 
locations  on  the  wing. 
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(c)  Deform  angle  is  the  local  deformation  angle. 

(d)  B  is  the  boundary  condition  (right-hand  side 
side  of  equation  57). 

2.  (a)  x  and  y  are  the  dimensional  location  of  the  control 

points. 

(b)  B  is  the  boundary  condition  (reference  equation  59) 

(c)  WW  is  a  dummy  variable  representing  part  of 
boundary  condition  term. 

(d)  BOD  DBLET  UPWSH  is  the  body  doublet  upwash. 

3.  (a)  Polynomial  coefficients  are  the  coefficients  which 

describe  AC  loading  distribution  on  the  wing 
(Reference  ^equation  54). 

4.  (a)  Pressure  distribution  for  wing  using  doublet 

paneling.  The  left-hand  column,  PSI,  are  the 
nondimen sional  chord-wise  locations;  the  top  row 
are  the  nondimensional  spanwise  locations  at  which 
the  pressure  is  computed. 

5.  (a)  The  next  output  is  the  sectional  coefficients 

for  the  deformed  wings  (reference  section  D,  above) . 

(b)  Next  is  the  total  coefficients  for  the  deformed 
wings  (reference  section  E,  above). 

G.  1.  XG  and  YG  are  the  x  and  y  locations  of  the  starting 

point  of  the  trailing  vortices  on  the  wing. 

2.  GMA  are  the  strengths  of  the  trailing  vortices. 

H.  Pressure  distribution  for  elevon  using  exact  functions 
(reference  section  4(a)). 

I.  This  section  contains  the  sectional  properties  for  the  elevon 
(reference  section  D) . 

J.  Total  force  coefficients  for  the  elevon  (reference  section  E) . 

K.  Conversion  to  Neilsen's  Data. 
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